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ABSTRACT 

We reexamine the impact of dark matter (DM) annihilation on the intergalactic medium, taking 
into account the clumping of DM particles. We find that energy injection from the annihilation of the 
thermal relic DM particles may significantly raise the gas temperature at high redshifts and leave a 
strong imprint on the cosmological 21-cm signal (6Tb > 10 mK), provided the particle mass is below 
~ 1 TeV. Further, we find that while the energy injection from DM annihilation could not alone 
complete the reionization of the Universe, it could make a significant contribution to the electron 
optical depth. 

Subject headings: cosmology: theory - early universe - galaxies: formation - galaxies: high redshift 
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1. INTRODUCTION 

In the standard cosmological paradigm, the non- 
baryonic dark matter (DM) particles make above 80% of 
the present matter density (|Spergel et al. 20071 ). How- 
ever, despite their high abundance, the DM particles 
have so far successfully resisted all efforts to detect 
them. While this non-detection has produced strong 
constraints on the nature of the DM particle and the 
strength of its interaction with the ordinary matter, 
there is still enough room for the non-gravitational 
processes involving D M to make a strong impact on the 
observable Univ erse (le.g. Berestrom fc Snellman 1988 1: 

Qin fc Wu 2"00lt iHennawi fc Ostriker 2 002: 

Chuzhov fc Nusser 20061: iBeacom fc Yuksel 200a 

Ascasibar 20071 : iMvers fc Nusser 20081 ). 

In this paper, we continue to explore the impact of 
the DM annihilation on the reionization history of the 
Universe. The products of the DM annihilation, which 
may include energetic photons and electron/positron 
pairs, are expected to couple to the ordinary matter, 
thereby raising the IGM temperature and the ionized 
fraction. In several previous papers, where this process 
was examined, it has been concluded that the effect of 
DM annihilation on the evolution of the intergalactic 
medium (IGM) during reio niza tion should be rather 
minor (iMapelh et al. 20061 : iFurlanetto et al. 20061: 
Ripamonti et al. 2007at Rip amonti et al. 2007b : 

Valdes et al. 20071) . The main reason for this con- 
clusion can be briefly described as following. Since 
the annihilation rate of non-relativistic DM particles 
is expected to depend primarily on their density, it 
should be the highest in the early universe, declining as 
(1 + z) 3 as the universe expands. Consequently, for the 
annihilation rate, which is low enough not to disturb 
recombination of Universe at z ~ 1000 violatin g the 
constraints set by the WMAP (jZhang et al. 2006[ L the 
energy input from DM annihilation during the epoch 
of reionization (z ~ 10) would be too low to make a 
significant impact. 

However, as we show in this paper, this conclusion may 
have to be changed, once the DM clumping is taken into 
account. While the DM density field is expected to be 
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nearly homogeneous prior and during the epoch of recom- 
bination, eventually the growth of density perturbation 
leads to formation of DM halos with constantly increas- 
ing density contrast. Consequently the DM annihilation 
rate declines much slower than the other processes that 
affect the IGM, such as inverse Compton (IC) scatterings 
and recombinations, and within the appropriate range of 
parameters DM annihilation can make a major contribu- 
tion to reionization, even if it were relatively unimportant 
earlier. 

In §2, we describe the evolution of the DM clumping. 
In §3 and 4, we estimate the impact of DM annihila- 
tion on the IGM and the cosmological 21-cm signal. We 
summarize our results in §5. 

2. DM CLUMPING 

The evolution of the DM distribution can be quali- 
tatively divided into several stages. In the early uni- 
verse (z <S 100) the DM density perturbations are gen- 
erally expected to be small and the clumping factor, 
Fcl = (Px) I '(px) 2 , close to unity. After a period of slow 
linear growth, density perturbations enter the non-linear 
regime, which is followed by formation of the DM halos. 
Unfortunately, the timing of this transition is uncertain. 
The amplitude of the primordial small-scale density per- 
turbations, which are first to enter the non-linear regime, 
cannot yet be probed observationally and may be esti- 
mated only by extrapolation from much higher scales. 
Further, the free-streaming of DM particles dilutes the 
perturbations by a factor that depends on a particle 
mass, which itself is an unknown quantity. Assuming 
the spectral index n — 1, which is favored by the infla- 
tion models, the non-linear regime should start around 
z ~ 50 - 70 (le.g. Green et al. 2005ft . However, one can 
not yet exclude the possibility that on small scales the 
primordial perturbation spectrum is drastically different 
from the assumed power law. If so, the non-linear regime 
may start shortly before reionization, or, at the other 
extreme, the small-scale perturbations may be already 
non-linear at the recombination epoch. 

After microhalos undergo the process of virialization, 
their internal density evolution slows down and their 
density contrast with respect to the average density of 
the expanding universe increases as ~ (1 + z)~ 3 . How- 
ever, when microhalos merge into larger objects, their 
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Fig. 1 — Thermal evolution of the IGM for f = 0, 0.1, 1 and 10 
(solid, dotted, dot-dashed and dashed curves, respectively). 

internal density is reduced by tidal forces and some 
microhalos may be completely evaporated. Numerical 
simulations indicate that after the merger the smaller 
halos typically lose most of their mass to the tidal 
forces of their neighbors or of the main halo. Nev- 
ertheless, the central regions of the microhalos remain 
relatively intact so that on average the tidal forces re- 
duce the dark matter clumping i s only by a factor of a 
few (le.g. Kazantzidis et al. 2004 iDiemand et al. 20061 : 
iGoerdt et al. 2007t iBerezinskv et al. 20071) . Later the 
formation of stellar disks may re sult in almost com- 
plete destruction of microhalos (|Angus fc Zhao 20071 : 
le.g. Zhao et al. 2005T) . but, since stars are expected to 
form in high numbers only at low redshifts, this is likely 
to happen only after the end of reionization. 

3. IGM EVOLUTION 
The energy input rate from DM annihilation is given 

by 

Ix = gfabsFd ( -f-^— ] (av)m x c 2 , (1) 

where f a bs is the fraction of the released energy absorbed 
by the IGM, (px) is the average DM density and mx is 
the mass of a single dark matter particle. It is gener- 
ally expected that the annihilation cross-section of non- 
relativistic DM particles scales inversely with their veloc- 
ity, so that their product (av) is a constant. The degen- 
eracy factor g equals 1, if the DM particle are identical 
(i.e., the particle and the anti-particles are the same) or 
1/2 otherwise. The value of f a bs depends on the model 
for the DM particle. Further, as shown by Ripamonti 
et al. (2007), it is likely to vary with redshift, possibly 
changing by more than an order of magnitude between 
z = 10 and 1000. However, to avoid restricting the cal- 
culations to a particular DM model, we shall have to 
assume that f a ts is a constant. 

We assume that, as is the case with most DM particle 
models, the energy from DM annihilation is injected into 
the baryonic matter primarily through collisions with 
high-energy electrons, which may be either the direct 
product of the DM particles annihilation or produced by 
secondary ionizations. A fraction of the absorbed energy 
is converted into heat, f Heat « (1 - (1 - x o.266)i.3i6^ 
where x is the ionized fraction, and the remainder is 



almost equally spli t between secondary ionizati ons and 
atomic excitations (IShull & van Steenberg 19851) . 

In the following calculations we assume that after en- 
tering the non-linear phase at z ~ 60, the DM clumping 
grows as F c i oc (1 + z)~ 3 . Prior to that epoch, we make 
a simplistic assumption that the density perturbations 
evolve as spherical top-hats. It should be noted though 
that the state of the IGM during reionization is almost 
insensitive to the assumptions we make for z > 60, since 
most of the energy from the DM annihilation would be 
released long after that epoch. 

The evolution of the IGM depends on the product of 
several uncertain parameters 

(UA ( (av) \ ( F cl (l + z) 3 \ / m x \~ X 
^ \ 0.1 J v 10- 26 cm- 3 S"V V 10 s J \GeVJ 

(2) 

To set a plausible range for the value of C, we have to 
consider separately the plausible values for f a f, s , (av) and 
F cl (l + z) 3 . For n = 1, F ci {l + z)~ 3 should of order 
5 • 10 7 (for comparison, F c i(l + z)~ 3 of the DM halo 
that forms at z — 30 is ~ 6 • 10 7 , assuming the NFW 
density profile and the concentration parameter c=4 that 
is characteristic of the newly formed halos). If, as usually 
assumed, the DM particles were in thermal equilibrium 
during their decoupling epoch, their annihilation cross- 
section is around (av) ~ 2 • 10~ 26 cm~ 3 s -1 . During 
the epoch of reionization f a b s is likely to be of order 
0.1 (jRipamonti et al. 20 07a). Thus we expect £ to be 
of order (mx /GeV) -1 . 

The WMAP observations s et an additional co nstraint 
on the DM annihilation rate (jZhang et al. 200"6h 

M< 0.8.10- 2 ^(^) cm" 8 "- 1 , (3) 

which for our choice of F c i translates into £ < 4. 

In addition to the energy injection from DM annihi- 
lation, the thermal evolution of the IGM is set by the 
adiabatic cooling, the change in the ionization state, the 
inverse Compton (IC) scattering of cosmic microwave 
background (CMB) photons and the radiative cooling 
processes 

dT 2f H _ 4T _ Tdx x Tcmb - T 

dt 3fc_Bn(l + x) 3t 1 + x l + x t 7 

( 4 ) 

where T is the IGM temperature, f 7 oc (1 + z) 4 is the 
IC cooling time-scale, fee is the Boltzmann constant, n is 
the baryon number density and Tcmb is the temperature 
of the CMB photons. The radiative cooling, L rac [, is 
dominated by collisional excitation of hydrogen atoms 
and hydrogen recombinations (we used the rates from 
Dalgarno & McCray (1972) and Seaton (1959)). The 
evolution of ionized fraction is dominated by radiative 
recombination and secondary ionizations induced by DM 
annihilation 

(rn\ 2 1 fionlx ,_\ 

_ = -a(T)*n+— , (5) 

where a(T) is the recombination rate, Eh — 13.6 eV 
is the hydrogen ionization potential and /j on ~ (1 — 
fHeat)/ '2 is the fraction of the absorbed energy that goes 
into ionizations. 

Figures Q] and [5] show the impact of the DM annihila- 
tion on thermal evolution and the ionization state of the 
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Fig. 2. — Same as Fig. [T]for ionization state. 
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Fig. 3. — The DM annihilation contribution to the electron 
optical depth. 

IGM for different values of (. Figure [3] shows the electron 
optical depth, Ar e , for the redshift interval 7 < z < 150, 
assuming the absence of other ionizing sources. When £ 
is close to its upper limit, we see that Ar e can exceed 
0.01. Further, if the DM clumping is a few times higher 
than we assumed, which is conceivable due to a large 
uncertainty in the small-scale DM distribution, then the 
Ar e can make a large fraction of the total electron op- 
tical depth, r e « 0.09 (jDunklev et al. 2008ft . However, 
we find that DM annihilation cannot alone complete the 
rcionization of the Universe. Since energy injection from 
the DM annihilation changes rather slowly, the end of 
rcionization (x = 1) must be preceded by a long period 
of significant partial ionization. Consequently, the en- 
ergy injection rate sufficient to complete rcionization by 
z ~ 7 results in the total electron optical depth being 
several times higher than allowed by the WMAP obser- 
vations. 

4. THE 21-CM SIGNATURE 

The next generation of radio-telescopes (LOFAR, 
MWA, SKA) should be capable of detecting the red- 
shifted 21-cm continuum from the epoch of reionization. 
Since the 21-cm signal is sensitive to the gas tempera- 
ture, ionized fraction and the UV intensity around the 
Lya resonance, DM annihilation, which affects each of 
these parameters, can produce a detectable imprint on 
the observed signal. To illustrate this, we calculate the 



evolution of the 21-cm spin and the differential bright- 
ness temperature, T 8 and 6Tb. 

The spin temperature is a weighted function of the 
CMB temperatur e, Tomb, a nd the kinetic temperature 
of hydrogen, T k , (jField 19581 ) 



r. 



Tcmb + VoTk + y c Tk 



1 



Va+Vc 



(6) 



where y c is the collisional coupling constant (see Zygel- 
man (2005) and Liszt (2001) for the contribution of neu- 
tral atoms and electrons to the collisional pumping) . The 
radiative coupling constant, y a , is 



Va = 



167r 2 !T,e 2 / 12 J 
27A w Tkm e c 



(7) 



where /12 = 0.4162 is the oscillator strength of the Lya 
transition, T* = hv*/k — 0.068 K, A w is the sponta- 
neous emission coefficient of the hyperfine transition, 
m e and e are electron mass and charge, and Jo is the 
intensity at Lya resonance, when the backreaction on 
the incide nt photons c aus ed by resonant scattering is ne- 
glected (iField 19581 : IChen fc Miralda-Escude 2004 ; 
Hirata 2006t~ IChuzhov fc Shapiro 2006 ; 

Furlanetto fc Pritchard 2006ft . For the unper- 

turbed IGM the backr eaction correction, S a , is 
(|Chuzhov fc Shapiro 2006ft 



S a = e -°-37(l+*) 1/2 T fc 



i + 2i 

T k 



(8) 



The Lya resonance photons can be produced by the con- 
ventional UV sources (i.e., stars and AGNs) or indirectly 
by DM annihilation. 2 

For unperturbed IGM the differential brightness tem- 
peratu re presently observ ed at the wavelength 21(1 + z) 
cm is (|Madau et al. 1997ft 



6T h « 0.03 K 



T s — T CMB 
T s 



(1-x) 



1 



10 



1/2 



(9) 



Figures [Hand [5] illustrate the evolution of T s and T^ in 
the absence of all other radiation sources. When the 
energy input from DM annihilation is relatively high 
(C <L l)i the evolution of T s can be distinctly different 
from other reionization scenarios. For instance, in mod- 
els where stars and miniquasars are the only radiation 
sources available, the decrease of T s due to the expan- 
sion of the Universe is expected to be reversed only a 
short time before reionization, at z < 20. By contrast, 
the DM annihilation may cause the decrease of T s be re- 
versed first already at z ~ 50, when the DM clumping 
enters the phase of fast non-linear growth. Subsequently 
T s can decline again until Lya photons are produced in 
large number either by stars or the line cooling in the 
IGM. 

Even a relatively low energy injection rate (C < 0.01) 
may produce an observable effect. In this case the DM 
annihilation does not create sufficient numbers of the 
Lya photons to decouple T s from Tqmb on its own. How- 
ever, it can still produce a noticeable effect on the gas 

2 The scatterings of non-thermal electrons, produced by DM 
annihilation, with hydrogen atoms, result in electron excitations, 
which i n ~ 5 out of 6 cases a re followed by the emission of the Lya 
photon llChuzhov et al. 20061) . 



4 



500 



200 



=,100 




10 15 20 30 50 70 100 150 



Fig. 4. — Same as Fig. [T]for the 21-cm spin temperature. The 
thick dashed line shows the CMB temperature. 
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Fig. 5. — Same as Fig. [T]for the 21-cm differential brightness 
temperature. 

kinetic temperature (Fig. [1} and, when the UV back- 
ground from other sources becomes sufficiently high, this 



would be reflected in the 21-cm signal. In this way, rais- 
ing the gas temperature by just a few degrees at z ~ 12, 
might change 5Tb by more than 100 mK, which is above 
the expected sensitivity (~ 10 mK) of LOFAR. 



5. DISCUSSION 

The high value of reionization optical depth, r e , mea- 
sured by WMAP, indicates that reionization might be- 
gin significantly earlier than the observations of the 
present-day galaxies would suggest, and that new types 
of radiation sour ces may be required to explain it 
(|e.g. Gnedin 2007t ). The high opacity of the IGM to the 
Lya photons at z > 6 further indicates that reionization 
proceeded over very extended period. Thus preionization 
produced by DM annihilation can fit well into the exist- 
ing constraints on the reionization scenarios. Naturally, 
not all of the DM candidates are suitable. The constraint 
on the relic abundance of the thermalized DM particles, 
(av) < 2 • 10~ 26 cm -3 s _1 , rules out most of the heavier 
DM candidates, such as the neutralino, whose annihila- 
tion rate would be too low. On the other hand, some 
of the lighter DM candid ates with mass below 100 MeV 
(|e.g. Hooper et al. 200l may fit the profile. 

DM particles, whose annihilation rate is too low to 
make a significant impact on the reionization, might still 
leave a potentially detectable signature in the cosmolog- 
ical 21-cm background, which is sensitive to the gas tem- 
perature. However, the heat input from DM annihilation 
may turn out to be secondary to other processes, such as 
X-ray absorption or scatterings o f UV resonance pho- 
tons fle.g. Venkatesan et al. 200ll lOh fc Haiman 20031 : 
IChuzhov fc Shapiro 2007t ICiardi k. Salvaterra 20071 ). In 
addition, given that during the reionization epoch the 
DM particles annihilation time-scale is expected to re- 
main nearly constant, the effect of the DM annihilation 
may be difficult to distinguish from the effect of DM de- 
cay. 
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